Electrodynamics of metallic photonic crystals and problem of left-handed materials 



A. L. Pokrovsky and A. L. Efros 
Department of Physics, University of Utah, Salt Lake City UT, 84112 USA 



O 



An analytical theory of low frequency electromagnetic waves in metallic photonic crystals with 
a small volume fraction of a metal is presented. The evidence for such waves has been obtained 
recently by experiments and computations. The cutoff frequency of these waves, ujq is studied. 
An analytical expression for the permittivity e is obtained and shown to be negative below wo- 
lf the crystal is embedded into a medium with a negative fj,, there are no propagating modes at 
any frequency. Thus, such a compound system is not a left handed material (LHM). The recent 
experimental results on the LHM are discussed. 
In his seminal work Veselago Q has shown that if in some frequency range both the permittivity e and permeability 
fi are negative, the electromagnetic waves (EMW's) propagate but they have some peculiar properties. All these 
■ properties come from the fact that vectors k, E, H form a left handed rather than a right handed set. It follows that 
the Poynting vector and the wave vector k have opposite directions. The materials with these properties are called 
C$ • the left handed materials (LHM's). 

The idea that a metallic photonic crystal (MPC) may be a technological base for the LHM ^|J|] appears as a 
result of the computational and experimental studies of a few groups Ji|-||] which have found the EMW's in the MPC 
propagating above a very low cutoff frequency. The MPC's they have considered are three- or two-dimensional lattices 
of thin straight metallic wires. Their discovery is very interesting and important because the waves propagate under 
' the condition fa/eoto 3> 1, where fa is an average conductivity, / being the volume fraction of a metal in the system. 
. This propagation must be due to the MPC structure, because it would be impossible in a homogeneous medium with 
the conductivity fa. 

Various groups obtained different cutoff frequencies and they understood them differently. The group of Soukoulis lq| 
qualitatively interpreted the effect of propagation in terms of waveguide modes, while the group of Pendry presented 
a completely original physical picture based upon a new longitudinal mode, called "plasma mode" . According to 
Pendry et al. M the resulting permittivity has a plasma-like behavior 
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however, the "plasma" frequency contains the light velocity and has a form 

§' 27rc2 r- 60 ^ (2) 

Q i <P\n{d/RY ~ ir&a ' 1 ' 

where d is the lattice constant, R is the radius of the metallic wires, a is the static conductivity of the metal. The 
same results for e and ui p have been later obtained theoretically by Sarychev and Shalaev ||. 

The San Diego group has accepted the "plasma model" and considered the negative e at ui < uj v as one of the two 
crucial conditions for creation of the LHM. To obtain negative fj, the split ring resonators (SRR's) are added to the 
MPC ||[lO|,[ll]] . The first observation of the negative refraction at the interface of this compound system and vacuum 
has been reported recently pQ| . The negative refraction is the most important manifestation of the LHM. 

We claim in this paper that the "plasma mode" in the MPC is in fact an EMW and that is why the experiments of 
the San Diego group cannot be simply understood using e of the MPC and fi of the SRR's. Our paper is organized as 
follows. First we discuss the arguments of Pendry et al. in favor of the plasma model and the theoretical approach 
by Smith et al. jq|. Then we derive and solve an exact dispersion equation for the cutoff frequency loq and find e(u>). 
The results are different from Eqs.(Q,^) but uiq is in a very good agreement with all computational and experimental 
data we are aware of. It follows that all groups discuss the same mode. The permittivity becomes negative at uj < luq. 
Then we show that the MPC does not support any waves if it is embedded in the medium with negative fi and discuss 
this result together with the experimental results of the San Diego group. 

The Pendry group has proposed that the plasma mode in the MPC appears as a result of an electromagnetic 
renormalization of the electron mass in the second order in 1/c 2 . Due to this renormalization the electron mass 
becomes ss 15 times larger than the mass of a proton. As a result, the plasma frequency of the metal shifts down into 
the GHz range. It is known, however, that the fields can be excluded from the interaction energy and the Lagrangian 
can be written as a function of instantaneous velocities and coordinates of the interacting charges keeping terms of the 
order of c~ 2 This so-called Darwin Lagrangian does not contain any mass renormalization due to the interaction. 
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It is also strange that one should care about both mass renormalization and plasma frequency in connection with a 
problem based upon Maxwell equations with a static conductivity. These equations contain neither electron mass nor 
plasma frequency, and the very concept of electrons is not important for them. 

The interpretation of the electric permittivity [Eqs.(]l],||)] as given by Smith et al. Q is in terms of conventional 
electrodynamics and it provides a reasonable basis for a discussion. As far as we understand, it is based upon the 
ansatz 

E ac = E + nrR 2 u:Lj, (3) 

where E ac is the electric field acting on a wire, E is the average field, j = aE ac is the current density, and L is the 
self- inductance of the wire per unit length. Using the ansatz (||) one can easily get Eqs.(Q,^). 

Now we derive an exact dispersion relation for the s-polarized EMW under the condition / <C 1 in the system of 
infinite parallel thin straight wires ordered in a square lattice. The electric field of the wave is along the wires (z-axis), 
while the wave vector k is in the x-y plane. Assume that the total current in each wire is Iq exp[— i(u)t — k • r,)], where 
Yi is the two-dimensional radius-vector of the wire in the x-y plane. The external solution for electric field E z of one 
wire with r,; = has a form 

E\ = (/ /W4)i4 2) (^/c), (4) 
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where Hq ' = Jq — iNg is the Hankel function which decays exponentially at Imw < 0. Neither Jo nor Nq has 
this important property. Here and below we omit the time dependent factor. The solution is written in cylindrical 
coordinates z,p,</> and it obeys the boundary condition = (i/uj)dE z /dp = Iopo/2irp at p = R. 
The electric field created by all wires is 

E z (r) = ^e^^- r) <\- c P^ (5) 

j 

where pj = y/Jx — Xj) 2 + (y — yj) 2 , summation is over all sites of the square lattice and the sum is a periodic function 
of r. 

The dispersion equation follows from the boundary condition |L3| that relates the total electric field at the surface 
of any wire / to the total current through this wire E z i = Iq exp (ik • ri)/a e firR 2 , where er e / = 2a J\(kR) / ' kRJq(kR) , 
k = (1 + i)/S, and 5 is the skin depth. At small frequencies, when S > R one gets a e f ~ <J- At high frequencies, when 
5 < 0.1R, one gets the Rayleigh formula a e f ps (1 + i)a5/R. 

Note, that the EMW exists mostly if the skin-effect in the wires is strong. Using Eqs.([l],^]) one can show that 
T/iOp — 5 2 /{R 2 In d/R), where S is the skin depth at uj — uj p (see also Ref. |TJ]). We show below that the exact 
solution has similar properties. That is why we mostly concentrate here on the case of the strong skin-effect. 

Finally, the dispersion equation for w(k) = (x — i"f)c/d has a form 

(x-i^e^^H^iz^) = (6) 



where z im = (x — i^l)\/l 2 + m 2 + (R/d) 2 , I and m are integer numbers, the small term (R/d) 2 under the square root 
is important only when I = m = 0. Taking real and imaginary parts of Eq.(^) one gets two equations for x and 7. 
In the continuum approximation one can substitute the summation by integration in Eq.(^) to get 

l-^+*^=0. (7) 

Equation (^) describes propagation of a plane wave through a homogeneous medium with the conductivity fcr e f, 
which is possible if fa e f/e^Lo <C 1. However, outside the continuum approximation there are propagating modes in 
the low frequency range fa e f/eoui 1. For these modes the fields are strongly modulated inside the lattice cell and 
E z (r) is close to zero near each wire so that the absorption is small. 

We begin with the frequency u>o which is the solution of Eq.(||) at |k| = 

( X -*7)£ ff o 2 W) = -i^- (8) 
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This mode is an eigenmode of the system and its frequency is the cutoff frequency for the EMW's. The numerical 
results for the real part of the frequency are shown in Fig. [l| One can see that x is °f the order of few units. The 
values of 7 are of the order of the right hand side of Eq.(||). Thus, x ^ 7 if /Ce/Ao^ ^ 1- 

In addition to the numerical solution we propose an approximation valid at very small /, when |ln/| 3> 1. We 
separate the term with I = m = and substitute the rest of the sum by the integral. Then 
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where we assume that 7 <C x- Here C is the Euler's constant. The second term at the left hand side represents the 
average field E z which can be found from the Maxwell equation by the following way. One can show that the average 
(or macroscopic) magnetic induction B is zero. Then V x B = 0, j z + £$dE z /dt — 0, and 



E z = 



iujeod 2 
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The second term in the square brackets of Eq.(||) is 4E z d/ hqIqc — —4d/(x ~ *7)- Thus, the expression in the square 
brackets describes deviation of the field acting on a wire from the average field and it is assumed that only one term 
of the sum in Eq.(^) makes this difference. This approximation is similar to the ansatz of Smith et al. The main 
difference between Eq.([3j) and Eq.(^|) which is crucial for the imaginary part of the frequency and also important 
for the real part, is the frequency dependence of a e f due to the skin effect. Say, in the work by Smith et al. Q 
5/R ps 7 • 10 -4 , so that the skin effect is very strong. 

Figure |l| compares our result for the real part of the cutoff frequency given by Eq.(||) with the results given by 
Eqs.(pph. In these calculations we assume that a e f is given by the Rayleigh formula so that the right hand side 
of Eqs.(|||)j) has a form a(l — i)y/x, where a = 2ceoR/ fdcrS and S is taken at u> = c/d. One can see from Fig.|l| 
that approximation Eq.(||) is much better than approximation Eq.(^). Both approximations coincide at small / 
and are accurate at extremely low values of / (~ 10~ 7 ) when the logarithmic term in Eqs.(|]j|) is very large. The 
computational and experimental data of Ref. Jg,|| are also shown at Fig. |l| and they are in a good agreement with 
Eq.(||). The results of Pendry group J?J (not shown) are also in a good agreement with Eq.(||). Thus, we can make the 
conclusion that the San Diego group, group of Pendry and Soukoulis group discuss the same mode but at different 
values of parameters and that our analytical theory describes the same mode as well. 

Now we find the component of electric permittivity e(u>) = e zz , which describes the s-polarized extraordinary waves 
in the uniaxial crystal. It is defined by the relation 
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where effective macroscopic conductivity a relates average current density j 2 to the average electric field E z by 
equation ~j z = bE z . To find a we introduce an external electric field E z e~ %ult . The average field E z is given by 
equation 



E z = £ z — 



d(x ~ ij) : 
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where the second term is the discussed above average field created by the wires. The boundary condition on a wire 
now has a form 



locjx ~ n)fJ-0 \ - rr(2) 
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Making use of Eqs.( |l2| , |l3| ) one can find a relation between the current and the average field which gives both a and 
e. Finally one gets 
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Assuming that the medium is transparent (7 — > 0, a — > 00) one can get the electric permittivity from the total energy 
density U of the electric field U = (l/2)E 2 d(ive)/du). We have checked that this method gives the same result for 
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e(u>). The expression in the square brackets of Eq.(p^) is the dispersion equation (||). One can see that Ree changes 
sign at uj = ujq and becomes negative at u> < loq, where ujq is the root of the dispersion equation (||). The derivative 
d(Ree)/du> at ui = ujq is the same as the derivative which follows from Eq.(jlj) at uj = ui p . However, far from the root 
of e equations (|lj) and (|j) differ substantially. 

To find w(k) one should solve Eq.(|^). For small |k| one can get analytical result uj 2 = c 2 k 2 + uj 2 ,, which is isotropic 
in the x-y plane. 

Now we discuss the possibility of creation of the LHM using the negative e of the MPC. Suppose that the wires are 
embedded into a medium with the negative magnetic permeability /i. One can see that in this case the propagation 
of any EMW is suppressed. Indeed, instead of Eq.(||) one gets equation 

STe^+^Ko (^J x 2 + y 2 + R 2 ) = | , (15) 



where cq = 1/ •v/j/Ifeo and Kq is the modified Bessel function. One can see that at |k| — all the terms on the left hand 
side of this equation are positive and real if Imw is small. Thus, if the right hand side is small, the equation cannot 
be satisfied. At small values of uj and |k| summation in the Eq.(|l5|) can be substituted by integration. Assuming that 
Reu) S> Imw one gets 

l + ^# = -^/Z. (16) 
uj z weo 

This equation does not have real solutions for uj(k). Thus, at negative /i there are no propagating modes at any 
frequency under the study. 

This result obviously follows from the fact that fieo < in the space between the wires. Therefore we get the 
plus sign in the left-hand side of Eq.([l6|) which forbids any EMW propagation. Thus, instead of the LHM we get a 
material without any propagating modes. 

Now we compare this result with the theoretical idea ||,|| to obtain the LHM, where negative e is created by the 
system of wires and negative [i is created in some other way. This idea is based upon the assumption that the negative 
e at uj < loo results from a "longitudinal plasma mode" . It is taken for granted that its frequency is independent of 
magnetic properties of the system, which is usually the case for plasmons. However, the mode discussed above is not 
a plasma mode (see also ) . One can show that this mode has zero average value of the magnetic induction B over 
the unit cell. In this sense this is indeed a longitudinal mode. But the average value of the magnetic energy, which 
is proportional to B 2 , is not zero and it is large. The physics of this mode is substantially related to the magnetic 
energy. That is why negative fi completely destroys this mode. It destroys also the region of negative e. In fact 
this could be predicted from the observation that ujq ~ c/d becomes imaginary at negative Say, one can see from 
Eqs . (0j2j) that at fx < one gets uj 2 < and e > at all frequencies assuming that T is small. 

Thus, we have shown that the simple explanation of the negative refraction in the compound system of the 
MPC and SRR's, based upon the permittivity e of the MPC and the negative permeability /i of the SRR's does not 
work because negative /i blocks propagation of EMW's in the MPC. The propagation observed by the San Diego 
group might be a manifestation of the remarkable conclusion of Landau and Lifshitz (See Ref. |l| p. 268) that (j,(uj) 
does not have physical meaning starting with some low frequency. Then, the explanation of the negative refraction 
in this particular system would be outside the simple Veselago scenario (see Ref. Q as an example). In this case, to 
explain the negative refraction one should use a microscopic equation similar to Eq. (^|) but with the SRR's included. 

Finally, an analytical theory of the low frequency EMW in the two-dimensional MPC is proposed. It is shown that 
the propagation of the low frequency waves is possible because electric and magnetic fields in the wave are strongly 
inhomogeneous inside the lattice cell and that electric field is small near the wires. If the dielectric part of the MPC 
has a negative fi, no waves can propagate through the system. We argue that the explanation of the experiment Ref. 
fiof is much deeper than it has been supposed before. 

We are grateful to A. Zakhidov and V. Vardeny who have attracted our attention to this interesting problem 
and to L. P. Pitaevskii and B. L. Spivak for important discussions. The work has been funded by the NSF grant 
DMR-0102964. 
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FIG. 1. The real part of the dimensionless cutoff frequency Xo as a function of the volume fraction of metal /. Solid, dashed, 
and dotted lines represent solutions of Eq.(|^), Eq.^, Eq.Q) respectively. The experimental data of Ref. [j^K*) are shown 
together with numerical data of Ref. Q (A, •, ■). On the main plot a — 0.024, d = 12.7fim for upper solid, upper dashed lines, 
and for the point A; a — 0.078, d — 1.27/im for middle solid, middle dashed lines, and for the point •; a = 0.246, d — 0.13^m 
for lower solid, lower dashed lines, and for the point ■. On the insert a = 3.4 ■ 10~ 4 , d = 8.0mm. 
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